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ABSTRACT 

The  principle  of  virtual  work  was  used  to  determine 
several  typical  difference  equations  governing  the  dis- 
placements of  a  taut  string  subjected  to  a  distributed 
transverse  load.   The  integrals  for  the  strain  energy  and 
the  external  work  were  approximated  by  both  the  trapezoidal 
rule  and  Simpson's  rule.   Further.,  the  integrand  in  the 
strain  energy  expression  was  replaced  by  two  different 
finite-difference  approximations.   For  comparison^  two 
typical  difference  equations  were  developed  from  the 
differential  equation  for  the  string  using  central  finite 
f^ffpT-pnfp.q  and  the  Hermitian  differencing  scheme.   The 
finite  difference  equations s    all  of  which  are  tridiagonal5 
were  evaluated  by  comparing  truncation  error. 
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I.   INTRODUCTION 

The  principle  of  virtual  work  is  a  very  useful  tool 
for  finding  the  displacements  of  an  elastic  body  produced 
by  given  forces  £l~]  .   In  this  thesis  the  principle  is 
used  to  determine  several  finite  difference  equations 
governing  the  displacements  of  a  taut  elastic  string  sub- 
jected to  a  distributed  transverse  load  per  unit  lengthy 
q.   The  strain  energy  of  the  string.,  U_,  is  given  by 


'L 


U 


dw 


2  Jo     \   dx 


dx 


(1.1) 


in  which  F  is  the  constant  internal  force  in  the  stringy 
w  is  the  displacement f    L  is  the  length,,  and  x  is  the  axial 
coordinate.   The  positive  directions  01    F >    w3    q5  and  x 
are  indicated  in  Figure  1. 


w 


Figure  1.   A  STRING  SUBJECTED  TO  A  TRANVERSE  LOAD 

AND  A  TENSILE  FORCE 
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The  principle  of  virtual  work,  is 

£U  =<fW  (1.2) 

where  fu  is  the  variation  of  U  with  respect  to  w  and  £w 
is  the  work  performed  by  q  moving  through  the  virtual 
displacement  <£w.   Hence 

<£W  =  J      q  <fwdx  (1-3) 

The  trapezoidal  rule  and  Simpson's  rule  are  used  to 
numerically  evaluate  the  continuous  integrals  in  Equations 
1.1  and  1.3  when  the  string  is  divided  into  k  intervals  of 
length  h.   The  displacements  of  the  string  are  defined  at 
the  end  points  or  nodes  of  these  intervals s    and  the 
continuous  first  derivative  in  U  is  replaced  by  a  finite 
difference  expression.   Several  approximations  to  the  first 
derviative  are  considered.   This  procedure  produces  a 
lengthy  equation  in  terms  or  trie  displacements  and  loads 
at  the  finite  number  of  nodes .   Taking  the  variation  of  this 
equation  with  respect  to  the  displacement  at  each  of  the 
nodes  3  in  accordance  with  the  principle  of  virtual  work., 
leads  to  a  set  of  simultaneous  linear  algebraic  equations. 
There  is  one  such  set  of  equations  for  each  numerical  scheme. 
Within  each  set  there  are  one  or  two  typical  difference 
equations.   The  schemes  are  evaluated  by  examining  the 
truncation  errors  associated  with  each  typical  equation. 


-*-The  variation  must  satisfy  the  geometric  boundary 
conditions _,  i.e..,  £w  =  o  at  x  =  o  and  x  =  L. 


A  second  procedure  for  finding  the  displacements  of 
an  elastic  body  is  to  use  a  finite  difference  approximation 
of  the  governing  differential  equation.   For  the  stringy 
the  differential  equation  is 

»  (§)  +  q  -  o  (1.4) 

Substituting  a  finite  difference  approximation  for  the 
second  derivative  leads  to  a  set  of  simultaneous  linear 
algebraic  equations  of  tridiagonal  form.   Two  typical 
finite  difference  equations  are  derived  and  compared  with 
the  typical  equations  developed  from  the  virtual  work 
principle . 
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II.   DEVELOPMENT  OF  THE  FINITE  DIFFERENCE  EQUATIONS 
USING  THE  PRINCIPLE  OF  VIRTUAL  WORK 

A.   THE  TRAPEZOIDAL  RULE  AND  SIMPSON'S  RULE 

There  are  many  formulas  available  for  numerically 
evaluating  integrals.   The  trapezoidal  rule  and  Simpson's 
rule }    two  of  the  most  commonly  used,,  were  chosen  for  this 
thesis.   Both  of  these  find  approximations  to  integrals 
by  passing  an  Interpolating  function  through  several  nodes 
and  integrating  the  interpolating  function.   The  trape- 
zoidal rule  uses  a  straight  line  a.s  the  interpolating 

function  between  x  and  x  .,,  where  x.  is  the  value  of 

s      s+15        s 

x  at  the   th  node,  and  is  given  as  J_  2  J 

/  ~  _  \ 

l  <£  .  J-  ! 

f(x)dx=2  f(xs+1)  -  f(xs)  -j2a^1 


in  which  x  <  T  <  x   .,  .   The  last  term  in  Equation  2.1 
s -  b "   s+1 

is  known  as  the  truncation  error.   The  integral  from  x=0 

to  x=L  is 

(2.2) 

f   f(x)dx  =  hJl-f(0)+f(x1)+...+f(xk_1)+if(L)l 

12    dx^ 
in  which  0   <     6  —  L 

In  Simpson's  rule  a  parabola  is  used  as  the  interpolat- 
ing function  passing  through  the  nodes  x   s   x  and  x    . 

S —  1     S         S  ~r_L 
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The  integral  from  x-0  to  x=L  is 
0 


J 


(2.3) 
f(x)dx  =  h  rf(0)+4f(x1)+2f(x2)+...+2f(xk_2) 


+4f(x    )+f(L)]  -   A  dViXl 
K_1       J     90   dx^ 

where  0  <  X  -  L  and  k  is  even. 


B.   DIFFERENCE  APPROXIMATIONS  TO  THE  FIRST  DERIVATIVE 

The  equation  containing  the  conventional  central 
finite  difference  approximation  to  the  first  derivative 
at  x   and  the  truncation  error  is  [_2J 

/dw\   =  1  (w    -  w   )  -  h2/dV)      (2.4) 

in  which  x   _  <"  Vs  <   x 

s-1-      s+1 

Consider  the  bcheme  in  which  the  displacements  are 
defined  at  the  nodes  and  the  first  derivatives  are  defined 
at  the  center  of  the  intervals  as  shown  in  Figure  2.   This 
scheme  is  referred  to  as  the  staggered  differencing  scheme 
or  half-station  scheme.   The  central  finite  difference 
equation  for  the  staggered  scheme  is  derived  in  Appendix 
B5  and  is 

i  (d  \ 

iS  dx~  .s+.2 


/dw\    =  1  (w    -  ws)  -  |£(l_w)  (2.5) 

I  dx/ _  ,  i   h  48  dx^  600.1 


where  x  t  T    ,ilx  ,,.   Hereafter,  unless  noted  otherwise 
s   °s+2   s+1  ' 

Xa   denotes  that  X   is  located  in  the  interval  x   ,  to  x 
s  s-1     s+1 

and  7*  04- i  is  in  the  interval  x  to  x    ,  and  Z  is  in  the 
S^2  s      s+1 

interval  x=0  to  x=L. 
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WjJ4, 


*S*i 


Figure  2.   THE  STAGGERED  SCHEME 

Another  finite  difference  equation  used  to  replace  the 
first  derivative  at  x  is  this  higher  order  expression 

/dw\  =  1  (w  in  -  w  J  +  h      (q    -q)+  7h%5^ 
Idx/    2h  v  s+1    s-ly   12F  VHs+l   *s-l'  3~60  (-^5) r 

(2.6) 
The  corresponding  higher  order  staggered  scheme  equation  is 

=  1    (w  ,  _  -  wj  +  _JL  (q„,n  -  q  )  +  7h  f  dW 

"2  u3  +  2 

(2-7) 

Equations  2.6  and  2.7  are  derived  in  Appendix  A. 


C.   DEVELOPMENT  OF  THE  TYPICAL  FINITE  DIFFERENCE  EQUATIONS 
The  first  typical  difference  equation  is  derived  using 
the  trapezoidal  rule  and  the  central  finite  difference  ex- 
pression.  Applying  Equations  2.2  and  2.4  to  Equation  1.1 
leads  to* 
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u  = 


Fh 


+p  K  -  ws-2)  -  £ 


Vdx 


^s-1 


+ 


1 

2h 


(w 


s+1 


-  W   -  )  - 


S-l 


h 

6 


6 


d3w 


"1  2 


dx 


3 


rc 


s+l  J 


+ 


dV 
,dx^, 


2 


+ 


? 
24 


2h 


d 


(ws+2~  ws) 
dw\2 

dxy^ 

2 


dx 


(2.8) 
Only  those  terms  in  which  w  appears  are  listed  in  Equation 
2.8  because  only  the  typical  equation  is  of  interest. 
Taking  the  variation  of  Equation  2.8  with  respect  to  the 


displacement  w  leads  to 
s 


lH   =  -  F  (wq  ?-2w  +w   -)-h2F  //^l 
<ws    4h    S~2    S   S+2    12  J  fe^ 


(2.9) 
/d3w) 


s+l 


Taking  the  variation  of  the  external  work  with  respect  to 


w  leads  to 
s 


&!      =     hq 


(2.10) 


Note  in  Equation  2.9  that  the  displacements  are  required 

at   n   and   „.   This  equation  does  not  take  into 
s-2     s+2         H 

consideration  the  fact  that  the  displacements  at   n  and 

^  s-l 

are  available.   Hence,  it  appears  reasonable  to  use 
s+l 

these  displacements  instead  of  those  at    and 

s-2    s+2 

Consequently,  let  the  interval  h  be  reduced  by  a  factor  of 
one  half  and  change  the  subscripts  accordingly  to  give  a 
set  of  equations  that  do  not  omit  any  nodes.   The  modified 
equations  are : 


14 


<^ws 


■  1   (w  , 

2h  v  s"l 

'd3w\ 


2w  +  w    ) 
s    s+ly 


dx^Rl 


(2.11) 


h 
2 


(2.12) 


It  can  be  shown  that  equations  2.11  and  2.12  are  identical 

to  the  equations  developed  using  the  trapezoidal  rule  and 

the  staggered  difference  scheme.   Using  the  staggered 

scheme  eliminates  the  necessity  to  halve  the  internal. 

As  an  example  of  the  derivation  of  the  higher  order 

equations,  a  typical  equation  is  derived  using  Simpson's 

rule  in  conjunctj.on  vn  th  the  hierher  order  staggered  scheme 

The  equation  for  U  based  on  Simpson's  rule  can  be  given  in 

the  form 

2       .  .  2 


U  =  ^   .  .  .  +  2  (• 

4/dw\  2 
h  L  d  [dxlZ 


S-p  S+2 


90 


dx 


J 

(2.13a) 


Equation  2.13a  could  also  have  been  given  as 


u  -25 


h 


.  .  .  +  4 


+  2 


/dw 
(  dx 


j     +  .  .  . 


2 

c 


dx 


(2.13b) 
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The  typical  equations  associated  with  Equations  2.13a  and 
2.13b  are  different.   The  equation  derived  from  Equation 
2.13a  is  referred  to  as  the  2-4-2  scheme  and  the  equation 
derived  from  Equation  2.13b  is  referred  to  as  the  4-2-4 
scheme. 

Substituting  Equation  2.7  into  Equation  2.13a  gives 


4     5 
7h     d  w ) 

5760  v  dx5^   1 
s  2  j 

+  24  (qs+l   <*s)  + 


h  (Ws  "  Vl)+2t  (%  -  Vl> 


+   4 


2 


4T   4  /dw\ 
(  dxJt 


h  L 
on 


—    (w   ,  -  w  ) 

h  v  s+1    s; 


d  w 


5760  (^5Jts+iJ 


+      .   .   . 


,w4 


/n  ili\ 


Taking  the  variation  of  U  with  respect  to  w  leads  to 

s 


e        =  "  —  (w    -  3w  +  2w  , _ )  -  $  (q   -  3q  +  2q    ) 
<U      3h  v  s-l   J  s     s+l;   36  VHS-1    s    s+ly 


+ 


4 
7h  F 

8640 


dJw 
.5 


(2-15) 


The  variation  of  the  external  work  with  respect  to  w  for 
the  integral  evaluated  by  Simpson  ' s  rule  with  the  2-4-2 
scheme  yields 

&  W    4h 


W 


3  ^ 


(2.16) 
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All  of  the  typical  difference  equations  that  were  derived 
using  the   principle  of  virtual  work  are  given  in  Tables  I 
and  II.   Table  I  contains  the  equations  derived  from  the 
trapezoidal  rule  and  Table  II  contains  the  equations  derived 
from  Simpson's  rule. 

The  truncation  errors  listed  in  Tables  I  and  II  are  only 
part  of  the  total  truncation  error.   These  truncation  errors 
are  associated  with  the  finite  difference  evaluation  of  the 
derivatives  in  the  strain  energy  integral.   Note  that  they 
consist  of  the  difference  of  two  derivatives.   As  a  con- 
sequence ^  they  can  be  converted  to  a  higher  order  of  h.   As 
an  example _,  from  the  first  equation  of  Table  I 

'*3~\  1  h2F     /dV\ 


hF    f~/d  w\  _   /d3w|  ~1  ^ 
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III.   THE  DIFFERENTIAL  EQUATION  APPROACH 

A.   CENTRAL  FINITE  DIFFERENCE  FORMULATION 

The  differential  equation  for  the  string  is  given  by 
Equation  1.4.   Using  the  conventional  central  finite  difference 
expression  to  evaluate  the  second  derivative  at  s  leads  to 

F(„s_i-2Ws+Ws+i)/h2  +  V  g  (±*}f  (3<1) 


B.   HERMITIAN  FINITE  DIFFERENCE  FORMULATION 

The  Hermitian  formula  for  the  second  derivative _,  as 
derived  in  Appendix  C,  is 


(w   -2w  +w    )  =  — 
V  s-1   s   s+i;   12 


*•     '   s-1  '  S       '  S+1J 


Applying  Equation  1.4  at     3      3    and     to  Equation  3.2 

s-1   s      s  +l 


gives 


F(w     -  2w   +  w    )  +  1  (q    +  lOq   4-  q    ) 

S-1       S      s  +  1     no  VHS-1         r      S+V 


12 


h^F  /dV) 

240  1^/5 


il,  =  0  (3.3) 


Equation  3.3  is  the  typical  difference  equation  for  the 
Hermitian  finite  difference  scheme. 
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Equations  3.1  and  3.3a  the  two  typical  difference 
equations,,  are  listed  in  Table  III. 
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IV.   CONCLUSIONS 

All  of  the  equations  derived  give  tridiagonal  matrices 
centered  about  the  principal  diagonal.   For  a  given  number 
of  nodes ,  a  computer  program  using  these  equations  should 
have  approximately  equal  storage  requirements,  computation 
time  and  round-off  error. 

There  are  two  ways  to  decrease  the  truncation  error. 
One  is  to  use  higher  order  difference  equations  and  the 
other  is  to  reduce  the  interval,  making  more  nodes.   Normally , 
both  of  these  cause  greater  round-off  error,  storage 
requirements ,  and  computation  time,  since  the  first  increases 
the  number  of  non-zero  terms  in  the  matrix  and  the  second 
increases  the  size  of  the  matrix. 

The  higher  order  equations  derived  for  this  thesis  do 
not  increase  the  size  of  the  displacement  matrix  but  do 
change  the  load  matrix  from  a  diagonal  to  a  tridiagonal 
matrix. 

The  object  of  any  numerical  analysis  method  is  to  obtain 
as  accurate  an  answer  as  possible  within  the  shortest  com- 
putation time.   Of  the  two  finite  difference  solutions  to 
the  differential  equation,  the  Hermitian  method  is  the  better 
of  the  two  from  the  standpoint  of  truncation  error. 

For  the  equations  derived  from  the  energy  solution, 
the  truncation  error  is  composed  of  three  parts.   Two  parts 
are  associated  with  the  evaluation  of  the  strain  energy  and 
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work  integrals  and  are  of  0(h)  for  the  trapezoidal  rule 

4 

and  0(h  )  for  Simpson's  rule.   The  third  is  the  one 

associated  with  the  evaluation  of  the  derivatives  in  the 
strain  energy  integral  and  is  listed  in  Tables  I  and  II. 

This  part  was  of  0(h  )  when  the  central  finite  difference 

4 

expression  was  used  and  of  0(h  )  when  the  higher  order 

finite  difference  expression  was  used. 

From  the  comparison  of  the  equations    made  in  this 
thesis  by  partial  truncation  error  it  would  appear  that  the 
energy  equations  derived  using  Simpson's  rule  and  the 
higher-order  finite  difference  equations  are  the  better 
equations.   However,  it  must  be  stressed  that  this  evaluation 
is  based  only  on  the  apparent  truncation  error  and  is, 
therefore,  not  conclusive. 

There  is  some  difficulty  in  establishing  the  total 
truncation  error  associated  with  the  equations  derived  from 
the  energy  approach.   Consequently,  it  is  recommended  that 
these  equations  be  evaluated  on  the  computer  so  that  a  com- 
plete comparison  can  be  made  not  only  of  truncation  error 
but  also  of  storage  requirements  and  computation  time. 
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APPENDIX  A 

THE  HIGHER  ORDER  FINITE  DIFFERENCE  EQUATION  FOR  THE  FIRST 

DERIVATIVE 


Writing  the  Taylor  series  for  w  _  and  taking  the 

s+1 

derivative  twice  gives 

(d2\A  =  /d2w  \    +   h/d3w\   +  h2/d  w  J   h3/  d^vA 

fev  B+1  ~u?va    Us;  s    2  i^vs  e  [  dx5/ 


+ 


(A.l) 


and   for  w      ,    gives 
s-1 


/ d  w  1         _  /d  w   j 
I  dx2  j  "I  d:,2  J 

»  '     S-]  R 


h/d3w  j       h2/d  w 
I  dx3/       2  I  d? 


m 


+ 


(A. 2) 


Subtracting  Equation  A. 2    from  Equation  A.l      gives 


,2 

d  w 


,,2 

d  w 


p  I  /       2 

dx    7  a+1   [dx     F  B-1 


A)  +h3/^w 


■  3 


^ 


3  /  s      3  I  dx- 


(A.3) 


Solving  this    f or/d  w  )  gives 

dx3 


d   yj 
dx' 


1 

2h 


d  w 


,2 

d  w 


2/  2 

dx  /  I  dx 

s+1 


s-1 


h2/A 
6     dx5 


+   . . . 


(A. 4) 
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Substituting  Equation  A. 4  into  the  central  finite 
difference  equation.  Equation  2.K,    for  the  non-staggered 
interval  gives 


/dw\    i 


1      12{.Ux2/s+1    U2;^      36o(dx5y 


(A.5) 


The  differential  equation  for  the  string  is 


d  w   — 


dx 


2 


F 


(A. 6) 


Substituting  Equation  A. 6  into  Equation  A.5  and  truncating 
gives 


(aw  » 
dx  J 


—  f  w  , ,  -w 
2h  *  s+1   s 


-q   J  +  7h  /d£w  ] 


12F 


(A. 7) 

Equation  A. 7  is  the  higher  order  finite  difference  equation 
for  the  non-staggered  scheme. 
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APPENDIX  B 

THE  DERIVATION  OF  THE  FORMULAS  AND  EQUATIONS  USED  WITH 
THE  STAGGERED  SCHEME 

A.   THE  CENTRAL  FINITE  DIFFERENCE  EQUATION  FOR  THE  FIRST 

DERIVATIVE 

The  Taylor  series  for  w    is 

s+1 

w  =  w      ,    +  h  l^i\  +  h2  (d  vj)        +     h3  (d\\ 

s+1         s+i       2  W)  s+l    H    ld^/s+i     4B    [to?/  s+L 


3«4 


/d\A  +  _h£     /d5w  )  + 

ld?fy        x        3840  \dl?/     - 

s+f  s+t 


(Bol) 


and  Taylor  series  for  w  is 

s 


/dw  \  +  h2  /Ar\  _   h3  /d3w\ 

(^js+i+H    ld#Js+i     W   {^5/  ^ 


,_2    /   2    \  ,.3 

s  s+i        ^ 


384 


UxV       x        3840  \dx^/ 


s+t  s+f  (B.2) 

Subtracting  Equation  B.2   from  Equation  B.l   gives 

'   w       -  w     =  h& )  +  h3  /d\  )         +  hL    /d^w)       +   .    . 

s+1       s  VWs+i       aP/s+i     1920l^/s+, 

(B.3) 
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/dw\ 
Solving  this   for  /  ~r~  1  gi 


lves 
s+i 


\dx/=    ,        h       s+1     s'        24  Ux3/  j      1920 1  dxV     , 

(B.4) 

By  truncating  after  the  term  in  w  the  central  finite 
difference  equation  for  the  first  derivative  for  the 
staggered  interval  is 

(  dx  J        l        h  (Ws+l   V  "  24  tdx3/?'  . 

where  x   -  T      ,  -   x  (B.5) 

s    vs+f    s+1  v    ' 

B.   THE  HIGHER  ORDER  FINITE  DIFFERENCE  EO.UATION  FOR  FIRST 

DERIVATIVE 

Writing  the  Taylor  series  for  w    and  taking  the 

s+l 

derivative   twice   gives 

/d2w\  /d2w)  h  /d  w  \  h2  /d  w\ 

\dx2)  =(dx7        1    MdxS/         +8-(a7V      x 


+  48 


s+i  s+f  s+f  s+i 

3  AM    +  ^to)     + 

U?/^    384  ( i^y  aai    •  •  • 

(B.6) 


s+f  s+^ 


and   for  w     have 
s 


2  °  r2 
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^8ldx5/  384    i?/ 

x  s+f 


+  .   .   . 


S+i 


(B.7) 


Subtracting  Equation  B.7   from  Equation  B.6   gives 


a2 
d  w 


s+l 


Solving  for 


(B.8) 


eives 


/a3„\ 


,2 
11    /d  w 


h        dx 


w  j         _     /d_w  J  _     h_  /d  w  )         + 

Vs+1      W/s  ^^/s+i 


(B.9) 

Substituting  Equation  Bc9   into  Equation  B.4   gives 


/dw\ 


=      i    (W  -W    )    -   ^ 

h        s+1      s7        24 


d  w 
dx' 


s+l     x         ' 


+ 


5760        dx 


d£w\  + 

dxV  s+i 


(B.10) 


Truncating  the   series    gives 


/dw\  =     1 

[axj       "    h 

s+^ 


(w    ,  _  -w    )    -  7- 
v    s+l      SJ      24 


/d  w  |  _  /  d   w 

K+1    fe 


+ 


7h 


4 


5760 


v  dxJ/rs+i 


(B.ll) 
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The  differential  equation  for  the  string  is 


jLs 

F 


m  ■■■ 

*    '    s 
Substituting  this  into  Equation  B.  11  gives 

(dw\  i  .      h      ,  N 

ta)        C     h    (Ws+rWs)+2^(qs+l-qS> 


4/5 
+  JOT 


5760  v  dxvr  1 

L   s+i  (B.12) 

Equation  B.12  is  the  higher  order  finite  difference 
equation  for  the  staggered  scheme. 
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APPENDIX  C 
THE  DERIVATION  OF  THE  HERMITIAN  FINITE  DIFFERENCE  EQUATION 

FOR  THE  SECOND  DERIVATIVE 


Adding  Equations  A.l  and  A. 2  gives 


4  /,6 

s 


+  h_   ill  +  .  .  .  (c.l) 


m 


Solving  Equation  C.l  for  /-t~*fI  Sives 

(aK)     if/d2w\    .  2  am  ,  AM  i-  hi  /d6wV 

Vdx-y  s      h-Max-/s.x        I  ^  /    3     [dx7s+j     ^2\dx/s 

(0.2) 

The  Taylor  series  derived  expression  for  the  second 

derivative  f  1 J  is 


dx2 


h^   b-1   s   s+l'    12  I  dx7   360  I  dx°/ 


(C3) 
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Substituting  Equation  C.2  into  Equation  C.3  gives 


AVj  .  i 


-  — o  (w   _  -2w  +w  ,  _  )  ~  — 
•  ^v  s-1   s   s+l;   12 


^  Ms-1 


+ 


h_ 
240 


^  A6 


rs; 


(C.4) 


Solving  Equation  C.4  for  (w   -2at  +V7  i   )  gives 

S  - 1     S    S~rl 


(w   n  -2w  +w  J  n  )/h2  =  -if  /^w 
^  s-1   s   s+l"     12   dx 


2 


+  10 


d  w 


s-1 


dx 


2 


hi  /A 

'240  I  7^ J  t. 


(C5) 


Equation  C„5  is  the  Herraitian.  equation  for  the  second 
derivative . 
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